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Abstract 

We consider a quantum system in contact with a heat bath consisting in an infinite 
chain of identical sub-systems at thermal equilibrium at inverse temperature /3. The time 
evolution is discrete and such that over each time step of duration r, the reference system 
is coupled to one new element of the chain only, by means of an interaction of strength A. 
We consider three asymptotic regimes of the parameters A and r for which the effective 
evolution of observables on the small system becomes continuous over suitable macroscopic 
time scales T and whose generator can be computed: the weak coupling limit regime 
A — > 0, T = 1, the regime r ^ 0, A^r and the critical case A^r = 1, r ^ 0. The 
first two regimes are perturbative in nature and the effective generators they determine is 
such that a non-trivial invariant sub-algebra of observables naturally emerges. The third 
asymptotic regime goes beyond the perturbative regime and provides an effective dynam- 
ics governed by a general Lindblad generator naturally constructed from the interaction 
Hamiltonian. Conversely, this result shows that one can attach to any Lindblad generator 
a repeated quantum interactions model whose asymptotic effective evolution is generated 
by this Lindblad operator. 



1 Introduction 

This paper is concerned with the study of the weak coupling limit, and variations thereof, 
of open quantum systems consisting in a small quantum system defined by a Hamiltonian 
/io on a Hilbert space coupled to a field or heat bath modelled by an infinite chain of 
identical independent n+l-level sub-systems on (i$n*TC, with n finite. The coupling between 
the distinguished system and the chain is provided by a discrete sequence of interactions 
of the small system with one individual sub-system of the chain, in the following way: if 
r > is a microscopic time scale, over a macroscopic time interval ]0, kr], k gN* , the small 
system is coupled with elements 1,2, ...,k of the chain, in sequence, for the same time r 
and with the same interaction of strength A. The interactions we consider are of the linear 
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minimal coupling type X]j=o ^* + cl'jj where the a*'s and a^'s are creation and 
annihilation operators relative to the levels of the sub-system and the I^'s are arbitrary 
operators on TCq. Such models of repeated quantum interactions are used in physics, e.g. 
in quantum optics, in the theory of quantum measurement or in decoherence. The lack of 
coupling, and thus of coherence, between the elements of the chain allows to expect that an 
effective continuous dissipative dynamics for pure states or observables on the small system 
of the form e*^ should emerge when the number k of discrete interactions goes to infinity 
and the coupling A with the chain elements is weak, in the familiar weak coupling regime. 
Recall that this corresponds to choosing t G M and considering N 3 k = t/X'^ so that the 
macroscopic time scale equals T = rt/A^. When r is fixed and A ^ 0, both k and T go 
to infinity as 1/A^. Moreover, in the setting adopted here, we have another parameter at 
hand which is the microscopic interaction time r of the small system with each individual 
element of the chain. It allows us to explore different asymptotic regimes, as r goes to zero 
as well, which characterizes the continuous limit, over suitable macroscopic time scales T. 

One goal of this paper is to establish the existence of effective continuous Markovian 
dynamics in weak and/or continuous limits defining three asymptotic regimes. We consider 
successively the effective Schrodinger evolution on the small system at zero temperature 
and, when the chain is at equilibrium at zero or positive temperature, the Heisenberg 
evolution of observables on the small system. While the existence of an effective dynamics 
obtained by a weak limit procedure (r = 1) is proven for a large class of time-independent 
Hamiltonian systems, as well as in certain time-dependent situations, see e.g. |Dlj . |D2j . 
|DSj . |LSj . |DJj . this question is not addressed in the literature for the case under study. 
Note also that a Hamiltonian formulation of our system necessarily involves a piecewise 
constant time-dependent generator. The analysis relies on the following property of the 
model, which is inherent to its definition. The effective dynamics on the small system of 
pure states or of observables from time to time kr is shown to be given by the k ^ 
power of a linear operator, which depends on the parameters A and r. This expresses the 
Markov property in a discrete setting. 

The first part of the paper is devoted to the usual weak limit regime A — > 0, r fixed 
and T = tr /}? — > oo, < i finite. We show the existence an effective dynamics driven by 
a r dependent generator which we determine. This first result is obtained by adapting the 
arguments developed in the study of the weak coupling regime for stationary Hamiltonians 
to our discrete quantum dynamics framework. The method is then extended to accommo- 
date the whole range r — > 0, AV ^ over macroscopic time scales T = t/{TX ) — > oo, 
which defines our second regime. This gives rise to an effective dynamics driven by a r 
independent generator we compute as well and to which we come back below. The analysis 
of these first two regimes is strongly related to regular perturbation theory in the param- 
eter A^r and we refer to these regimes as perturbative regimes. Technically, the study of 
the second regime relies on an asymptotic analysis in the two parameters A and r of the 
discrete evolution of our system. The divergence of the macroscopic time scale imposes, as 
usual, some renormalization of the dynamics by the restriction of the uncoupled dynamics. 
Finally, note that in the second regime, the interaction strength A is not required to go 
to zero and can even diverge. The common feature of the generators of the dynamics of 
observables obtained in these two regimes is that they commute with the generator i[/io, •] 
of the uncoupled unitary evolution restricted to Ti.Q. In other words, the corresponding 
effective dynamics admits the commutant of /iq as a non trivial invariant sub-algebra of 
observables. This property is well known in the weak coupling regime for time-independent 
Hamiltonians, [m], |EH1, |DJj. 



2 



Our primary motivation actually comes from the recent paper |APj where such repeated 
interactions models are shown to converge in some subtle limiting procedure to open quan- 
tum systems with a heat bath consisting in continuous fields of quantum noises, at zero 
temperature. These limiting systems give rise in a natural and spontaneous way to effec- 
tive dynamics on the Hilbert space TCq of the small system governed by quantum Langevin 
equations. The above mentioned limit involves at the same time the time scale r, the 
strength of interaction A as well as a notion of spacing between the sub-systems forming 
the chain in an intricate way. While reminiscent of weak coupling methods in spirit, the 
limiting procedure of |APj is nevertheless distinct from the weak coupling limit. Indeed, 
while r — > 0, the product A^r is kept constant in |APj . which leads us beyond the perturba- 
tive regime. Another goal of the present work is to consider as our third regime the critical 
scaling A^r = 1 in our repeated quantum interactions model and to derive an effective 
evolution of observables for a chain at inverse temperature /?. The relevant macroscopic 
time scale in this regime is T = t/(A^r) = t, which is finite. 

With this scaling, we show that an effective Heisenberg dynamics for observables on 
Tio emerges at any temperature. It is generated by a general Lindblad operator whose 
dissipative part is explicitly constructed in terms of the V^ 's defining the coupling in the 
Hamiltonian, whereas its conservative part is simply i[/io,-]. At zero temperature, we 
recover the effective Heisenberg dynamics of observables on Tio of |APj obtained by means 
of quantum noises. At positive temperature, our generator coincides with a construction 
proposed in |LMj for certain models using an a priori modelization of the heat bath by some 
thermal quantum noises, generalizing those used at zero temperature. For any temperature, 
the effective dynamics is distinct from that obtained in the previous two perturbative 
regimes. In particular, the generator obtained does not commute with i[ho, ■] anymore, 
the generator of the uncoupled evolution restricted to Tio. Hence, there is no obvious sub- 
algebra of observables left invariant by the effective dynamics of observables. The analysis 
of this critical case makes use of Chernoff 's Theorem, rather than perturbative methods. 

Let us compare the generator of the effective dynamics of observables obtained in the 
regime r — > 0, tA^ — > 0, and the general Lindblad operator obtained as r — > with tA^ = 1. 
In the former case, the generator is obtained from the dissipative part of the Lindblad 
operator of the latter case by retaining its diagonal terms only with respect to the spectral 
decomposition of the uncoupled evolution restricted to TCq. Or, in an equivalent way, by 
performing a time average of the Lindblad operator with respect to the uncoupled evolution 
restricted to TCq. This defines the so-called ^ operation that makes the commutant of /iq 
invariant under the effective dynamics in the regime r ^ 0, rA^ (and in the weak 
coupling regime as well). Our results show that the # operation is present as long as 
A^T 0, whereas it disappears in the critical regime rA^. In other words, in the regime 
T — > 0, tA^ ^ 0, a non-trivial distinguished invariant sub-algebra of observables exists, 
whereas in the critical case r — > 0, rA^ = 1, there is a priori no sub-algebra left invariant 
by the effective dynamics, since its generator takes the form of a generic Lindblad operator. 

We finally note here that from a practical point of view, the modelization of the dy- 
namics of observables (or states) of a small system in contact with a reservoir at a certain 
temperature often starts with a choice of a certain Lindblad generator suited to the physical 
phenomena to be discussed. Our analysis allows to assign to any Lindblad generator a sim- 
ple model of repeated quantum interactions, with explicit couplings constructed from the 
Lindblad generator, whose effective dynamics in the limit r — > 0, A = l/y/r, is generated 
by the chosen Lindblad operator. 
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The paper is organized as follows. The general setup and definition of the model are 
provided in the next section. Section 3 is devoted to the analysis of the weak limit of the 
model at zero temperature, in the Schrodinger picture. Our main results in this setup are 
expressed as Corollarv 13.21 for the weak coupling regime and Corollarv 13.31 for the regime 
A^T — >■ 0, T — > 0. This section also contains the technical basis underlying our perturbative 
analysese in both the Schrodinger and Heisenberg pictures. The main technical result, 
of independent interest, is actually valid in a Banach space framework and is stated as 
Theorem 13.11 The positive temperature case, in the Heisenberg picture is dealt with in 
Section 4. The generators of the effective dynamics of observables in the two perturbative 
regimes are given in Theorem 14.11 and Corollary 14.11 The analysis of the critical regime 
A^T = 1 is presented in Section 5, for both the Schrodinger and Heisenberg pictures. Section 
6 is devoted to a thorough analysis of the first non-trivial case where both the small system 
and the elements of the chain consist in two- level systems. 



2 A Repeated Interaction Model 

Consider the following setup to start with. Our small system, described by the Hilbert 
space 7^0 of dimension d + 1 > 1 and a self-adjoint Hamiltonian /iq, interacts with an 
infinite chain of identical finite dimensional sub-systems modelling a field or heat bath, 
by means of a time dependent Hamiltonian. The total Hilbert space is TCq (8> Ti., where 
n = 0j>iC"+\ n > 1. We win call the jth Hilbert space C^+^ = C"+\ the Hilbert space 
at site j, j = 1,2, ■ ■ ■ and, following the usage when n = 1, we will call the subsystem at site 
j the spin at site j. We adopt the following convenient notations used in jAPj . The vacuum 
Q H is defined as the infinite tensor product of the vacuum vector u; = ( • • ■ l)^ 
in C"+\ 

n = uj^LO^uJ^ ■■■ G C"+i C"+^ C"+i • • • . (2.1) 

Denoting the ith excited vectors Xj = (0 ••• 1 ••• 0)"'", where the 1 sits at the 
ith line, starting from the bottom, i = 1,2, - ■ ■ ,d, the corresponding excited state at site 
j > 1 is given by 

Xi{j) = a; (gi • • • a; (gi Xj (gi u; • • • , (2.2) 

where Xi sits at site j > 1. More generally, given a finite set 

S = {{ki,ii),{k2,i2), - ■ ■ ,{km,im)} C (N* X {1, 2, • • • , d})™ with all fc/s distinct, (2.3) 

we define Xs as the vector given by an infinite tensor product as above, with ijth excited 
vectors Xi.{kj) at all sites kj > 1, j = l,---,m, and ground state vectors oj everywhere 
else. This construction together with the vacuum Q. = Xijf yield an orthonormal basis of 
TL, when S runs over all finite sets of the type above. 

Let us introduce creation and annihilation operators associated with the vectors Xi{j). 
Let ai and a*, i = l,2,---,n, denote the operators corresponding to {u!,xi, - ■ ■ ,Xn} in 
C"+\ i.e. such that 

aiXi = uj, Qiio = GiXj = 0, if j 7^ i, 

a*oj = Xi, a*Xj=0 for any j = l,2,---,n. (2.4) 



4 



Note that these operators do not coincide with the famihar creation and annihihation, 
however, for i fixed, they satisfy the anti-commutation rules when restricted to the two di- 
mensional subspace < uj,Xi > and are zero on the orthogonal complement of this subspace. 
Then, for j > 1, the operators ai{j) and ai{j)* on 7i are defined as acting as and a* on 
the jth copy of C"+^ at site j, and as the identity everywhere else. Therefore, when acting 
on different copies of C""*"^, these operators commute. In keeping with the notations for 
the reservoir, we introduce a basis of eigenvectors of for Hq of the form 

{u,xi,X2-,- ■ ■ iXd}, where d = dimiT-Lo) — \. (2-5) 

Note that d ^ n in general, but we shall nevertheless use sometimes the notation a;(0) and 
{xi{0)}i=i^2, --,d to denote these vectors. No confusion should arise with vectors of H above, 
since we labelled the sites of the spins by positive integers. In some cases, Hq will be an 
infinite dimensional separable Hilbert space, which corresponds formally to d = oo. 

Our formal time dependent Hamiltonian H(t, A) onHo'S>H has the form 

H{t, X)=Ho + Hf + XHi{t), (2.6) 

where 

oo n 

Ho = ho^I, HF = Y,^^^^MjTai{j), with (5^ G M, (2.7) 

j=i i=i 

and, for t G [r(A; — 1), Tk[, 

n 

Hi{t) = Y^V*(^ ai{k) + Vi(^ ai{k)* = I{k), (2.8) 

i=l 

where the VJ's and ho are bounded operators on Ho, in case Ho is a separable infinite 
dimensional Hilbert space. These operators describe the interaction between the small 
system with the different levels of energy 5i of the spin at site k, during the time interval 
]r(/c — 1), rk] of length r. The form of Hp makes it an unbounded operator, but, as we will 
see in the sequel, we will only make use of the unitary evolution it generates and, moreover, 
it will always be sufficient to work with subspaces containing finitely many excited states 
only 

In order to make the notations more compact, we introduce vectors with operator valued 
entries that allow to get rid of the indices i = 1, - ■ ■ ,n. Let 

a(j)» = (ai(i)« a2(j)» ••• a„(j)» )^ (2.9) 
V^ = {VI Vl ■■■ Vl) (2.10) 

where " denotes either nothing or *. Then, using the rules of matrix composition, we can 
write 

n 

F«i ® a(i)fc = ® (2-11) 

i=l 

so that we can rewrite the interaction Hamiltonian for t g]t(A; — l),rA;] as 

I{k) = V* ®a{k) + V ®a{k)*. (2.12) 
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Similarly, with 

a{j)^a{j) = {a,{j)ia,{j) a2{j)Wj) ••• «n(j)«a„(j) f (2.13) 
5 = {5i 62 ••• 6n), (2.14) 

we can write 

HF = l0Y.6a{jra{j). (2.15) 
i>i 

We will denote the corresponding evolution operator between the time T{k — 1) and rk 
by Uk, so that 

Uk = e--(^o+/^F+A/(fc))^ (2.16) 
and the evolution from to rn is given by 

U{n,0) = UnUn-i---Uk---Ui. (2.17) 
Although not explicited in the notation, the operator U{n, 0) depends on A and r. 

We will first be interested in the weak coupling limit of this evolution operator charac- 
terized by the familiar scaling 

n = t/\^ , A ^ and r fixed. (2.18) 

Hence, the macroscopic time scale T is given by 

r = rn = Tt/A^ ^ 00. (2.19) 

Note, however, that in contrast with the usual set up, we have here a non-smooth time 
dependent Hamiltonian H{t, A). 

Remarks: 

i) In order not to bury the main points of our analysis under technical subtleties, we have 
chosen to work in a simple framework where all relevant operators are bounded or matrix 
valued. Nevertheless, some of our results below hold if we consider our heat bath to live 
in a tensor product of infinite dimensional separable Hilbert spaces and make further as- 
sumptions so that the field Hamiltonian and interaction are bounded. 

ii) In some cases we shall allow TIq to be a separable Hilbert space. This will be explicitly 
stated in the hypotheses. Otherwise, we will work on the model defined above, under the 
general assumption 

HO: The Hamiltonian is defined on the Hilbert space ® where = C*^^^, 
Ti = ®j>iC"'*'^, for d, n finite, and is given by (|2.6|) . (|2.7|) . (|2.8j) . The evolution it generates 
is given by ((TT7|) and (|TT^ . 
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3 Weak limit of the Schrodinger representation at 
zero temperature 

As a warm up, and in order to derive some preliminary estimates, we prove here the exis- 
tence of the weak hmit for our model at zero temperature in the Schrodinger picture, and 
compute this limit. We first prove a key lemma that reduces the computation of the pro- 
jected part of the evolution U{n,0) l\2.n\i to the n power of a single matrix. Then we 
perform a general analysis of large powers of operators based on perturbative expansions 
which appear in the computations of weak limits. These technical results are expressed in 
Proposition 13.21 and Theorem 13. II under different sets of hypotheses. Their applications to 
our model are given in Corollaries 13.21 and 13.31 



3.1 Markov Properties 

Let P be the projection from 7io 7i to the subspace Tio CQ defined by 

P = I®|Q)(Q|. (3.1) 

The object of interest to us in this Section will thus be the limit 

limP(C/(^/A^O))P, (3.2) 

as an operator from P Hq (dH to P TCq Ti, identified with Hq, the Hilbert space of the 
small system. 
Note that 

Uj = e-'^^^e-'^^\ (3.3) 

with 

Hj = ho®l + l® Sa{j)*a{j) + X{V* a{j) + V(g, a{j)*) 

Hj = I(^"^6a{k)*a{k), (3.4) 

two operators that commute. 

We observe the following property of products of operators Uk, which shows the Marko- 
vian nature of the reduced evolution. 

Lemma 3.1 Let us write the restriction ofUj to 'Ho'S>C^^^ as a block matrix with respect 
to the ordered basis of Tio ® C"^"*^ 

{lO (Si UJ, Xl LO, ■ ■ ■ ,Xd'Si UJ, 

uj (E) xi,xi (g) xi, - ■ ■ Xd(E) xi, 

(3.5) 

UJ (g) Xn, Xi® Xn, - ■ ■ Xd® X„} 

as 

where A is a (d + 1) x (d + 1) matrix, B is [d + 1) x n{d +1), C is {d+ l)n x (d + 1) and 
D is n[d+ 1) X n{d+ 1). Then, for any m > 0, 

PU{m, 0)P = A"'0 \n){n\ ~ A"". (3.7) 



7 



Proof: Follows from the fact that 

Uj{i(^\n){n\) = e-'^^^{i^\n){n\) (3.8) 

where, HHq B v = ^0^(0) + X^iLi 'i^iXi{0) ^ v, 

n+l 

e-""^v ®n = Av®Vt + Y^{Cv)^ (g) Xi{j), (3.9) 

i=l 

where {'w)i denotes the i'th component of the vector w. Hence, due to the fact that different 
C//s act on different C^+^'s, 

m{d+l) 

UraUm-l---UiV®Q. = A^V®Q.+ ^ Wi®Xs,, (3.10) 

i=l 

where Wi are some vector in C^^^ and the excited sets Si are never empty. Therefore their 
contribution vanishes in the computation 

(1 \VL){VL\)UmU^-i ■■■Uiv^n = A'^v^ (3.11) 



As is easy to check along the same lines, in case Tio is infinite dimensional, we can 
generalize the above Lemma as follows. 

Lemma 3.2 Let TCq be a separable Hilbert space and /iq, Vj, j = 1, • • • ,n be bounded on 
Ho- We set 

Pj = \xj){xj\ : C"+^ ^ C"+\ j = 1, • • • ,n, Po = \uj){uj\ and Qo = I- Pq, (3.12) 
so that 

no C"+i = (Tio PoC"+^) e {Tio QoC"+^) ~ {Tio C) © (Ho ® C"). (3.13) 
We can decompose 



where A : Hq ^ Hq, B : Ho ^ ^ Ho, C : Hq ^ Hq and D : Hq ^ Ho C". 
Then, for any m > 1, 

P?7(m,0)P = ^™® (3.15) 

The above Lemmas thus lead us to consider a reduced problem on Tio- We need to 
compute the matrix A in the decomposition (|3.6|) of g-^'^i^^o+Sa*a+x{v*a+Va*)) ^ where we 
dropped the indices j, the I and the (8) symbol in the notation. Recall however that a 
summation over the excited states of Hp is implicit in the notation. 
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3.2 Preliminary Estimates 

In order to apply perturbation theory as A — > and, later on, in other regimes involving 
T ^ as well, we derive below estimates to be used throughout the paper. 
We rewrite the generator as 

H{X) = H{0) + XW, with H {0) = ho + 6a* a and W = V*a + Va*. (3.16) 
With a slight abuse of notations, the projector P takes the form 

P = I-a*a. (3.17) 
We can slightly generalize the setup and work under the following hypothesis: 

HI: 

Let P be a projector on a Banach space B and H{X) be an operator in of the form 

H{X) = H{0) + XW, (3.18) 

where H{0) and W are bounded and < A < Aq for some Aq > 0. Further assume that 

[P, H{0)] =0 and W = PWQ + QWP where Q = I-P. (3.19) 

We consider 

Ur{X) = e-'^^^^\ (3.20) 

For later purposes, we also take care of the dependence in r of the error terms. As this 
parameter will eventually tend to zero in some applications to come below, we consider the 
error terms as both A and r tend to zero, independently of each other. We have a first easy 
perturbative result 

Lemma 3.3 Let HI be true. Then, as X and r go to zero, 

p^-ir{H{G)+\w) p ^ p^-iTH{o)p ^ x^PG{t)P + PO(A^r^)P, (3.22) 

where 



F{t) = XI H{Q)'^^WH{QY 



n>l m,SFi 



mi-\-7n2=Ti—l 



r dsie''^^^°^We-''^^^^^ (3.23) 
Jo 



g{t) = Y.^^^ Yl H{or'WH{or^wH{or^ 

n>2 ^' mjen 



Jo Jo 



Moreover 



-^G{t) = -iH{0)G{T) - iWF{T), G(0) = (3.25) 

P(_r) = -e^^^(°)p(r)e^^^(°) (3.26) 

G(-r) = -e'rHiO)Q^^yrH(0) ^ H (0) p ^^yr H (0) p ^^yr H (0) _ (3 37) 
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Remark: Formula H3.21() is true without assuming that W is off-diagonal with respect to 
P and Q. 

Proof: First note that Ur{X) = e~'''^^^^^ is analytic in both variables A and r in C^. Then, 
we compute the exponential of —it times H{X) as a convergent series. Consider terms of 
the form 

n-l 



{H{0) + XW)'' = HiO)"^ + X^H{0)''WH{0)''-^-^ 

+ Hior^WHior^WHior--' +0{X^C''). (3.28) 



The error term in C" comes from the boundedness of the operators involved. Multiplication 
by (— ir)"/n! and summation over n > yields the first result with our definition of F{t) 
and G{t). The second result follows from taking into account that PWP = QWQ = 0, 
hence only the terms with an even number of W^s survive and we get 

p(jm + mip^ (3.29) 

n! 



P 



n! ^ n! \ / ■) 



P. 



T=0 



The overall error in r^A^ comes from the fact that it takes at least four terms in (|3.28j) to 
get a contribution of order A"^. The computation above was conducted to order A^ because 
of the scaling 1)2. 18() . The order A term F{t) doesn't contribute, being off diagonal with 
respect to P. 

An alternative derivation of a perturbation series of e-*^(-f^(0)+^^) in A yields the other 
expressions for F{t) and G{t). It is obtained via Dyson series in the familiar interaction 
picture. We have the identity 

d_ 

Introducing 

e(A,T) = e-^^(0)e— (/^(0)W)^ (3,31) 
this operator satisfies 

z^e(A, r) = Ae*"-^W VFe-^"^(°)e(A, r), G(A, t)\^^^ = I. (3.32) 
Hence we have the convergent expansion 

e(A,T) = V(-iA)" / dsi ds2--- dsne'^'^^^^W X 

Xe~*(''i~'*2)J^(0)|yg-i{s2-S3)J?{0) . . . ^-is„-i-s„)H(0)^r^-iSnH{0) ^ (3.33) 

Therefore, focusing on the terms of order A and A^, we get the alternative expressions for 
F{t) and G(t). 
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The differential equation yielding G{t) as a function of F{t) follows from explicit com- 
putations on the expressions above, as the identities for r i-^ — r. I 

Let us give some more properties of the expansion of Ur{X) for A > small, r > in 
the Hilbert space context that will be used later on. 

Corollary 3.1 Assume B is a Hilbert space, H(0), W and P are self-adjoint and A,r are 
real. As X ^ 0, the operator Ur{X) = e^*'^^^'^) satisfies 

Ur{\) = e-*^-^(°) + AF(r) + X'^G{t) + 0{X\^) (3.34) 

Ur{X)'^ = UriXr = U.riX) 

= e^^^(°) + AF(-r) + X'^G{-t) + 0{X\^) (3.35) 

with the identities for all t GM. 

F{-t) = F*{t) (3.36) 
G{-t) = G*{t). (3.37) 

Proof: Follows from the fact that H(X) is self-adjoint. I 



3.3 Weak Limit Results 

The technical basis underlying all our weak limit results is contained in the next two 
Lemmas and the Proposition following them. They are stated in a general framework that 
will suit both our analysese of the Schrodinger and Heisenberg representations. This is 
why we use independ notations. 

Lemma 3.4 Let V{x), x E [0,Xo), and R be bounded linear operators on a Banach space 
B such that, in the operator norm, V{x) = F(0) + xR + 0{x'^), and V{0) is an isometry 
which admits the following spectral decomposition 

r 

V{0) = '^e~'^^Pj where r < oo, Ej e R, {e"^-^^}j=o,-,r distinct. (3.38) 
i=o 



Let h = YJ'j=QEjPj so that 1/(0) = e"*'* and 



if j 



J=Y,»jkPjRPk where a^k = \ -f''' _ (3.39) 



j,k=0 



ie^^ if j = k. 



Then, for any < t < to, where to finite, and t/x E N, 

||y(x)^ -e-^('^+^'^)^|| =0(x), as x ^ 0, s.t. t/x E N. (3.40) 
Remarks: 

i) Expressing the projectors Pj by Von Neumann's ergodic theorem as 

p^ = ;^^^NY.^^''''^^^)r (3-41) 

n=0 

shows that they are of norm one. 

ii) The operator J = J{R, h) is defined as the solution to the equation ()3.43|) . This equation 
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is a particular case of i e'^^^Xe ds = Y which is solved in a similar fashion. 
Proof: With m = t/x G N, 

m—l 

V{xr - e-^('^+^^)"' = V{x)''{V{x) - ^-i(h+xJ)y-i{h+xtj)^-^-k ^3^42) 

fc=0 

where, by hypothesis and Lemma 111, SI 

V{x) - e-*('^+^-^) =x(^R + ie-'^' e'^'Je-'''' ds^ + 0{x'^). (3.43) 
Moreover, note also 

\\V{x)\\ = l + 0{x), ||e-^('^+^^)|| = l + 0(x). (3.44) 

Our definition (|3.39j) of J is designed to make the term of order x in (|3.43j) vanish. There- 
fore, there exists positive constants cq, ci such that we can estimate for any <t < to < oo 

m—l 

\\V{xr - e-^(''+^-^)™|| < cx^ ll^(^)f l|e"^^''+'^-^^ir"^"'' (3.45) 

k=0 

< cox^m{l + cox)™ < cotxe^ H^+cox) < xcotoe"'^'' = 0{x). ■ 

It will be necessary to control the dependence of such estimates on a parameter r ^ 
later on. This will cause no serious difficulty, since all steps are explicited in the argument. 
To achieve sufficient control in r, we need to revisit the proof of a well known lemma, which 
holds under weaker hypothesese than ours, see Davies |D2j . 



Lemma 3.5 Let e^*^ = Yl^j=o^~^^^ isometry it A'. !^8\} on the Banach space B and 

let K he a hounded operator on B. There exists a constant c depending on r and to only, 
such that for any t £ [0, to]? ^0 finite, 

_ ^-...|| < ^^Wmi + IlKlDe^ll^ll^o ^^^^ ^3^^^^ 

mfj^fc \ Ej — Ek\ 

where K* = YJ'j=o PjKPj = limy^oo ^ lo e'^'^Ke-^'^ ds. 

Remark: The expression of K"^ as a Cesaro mean is a classical computation which shows 
that \\K*\\ < \\K\\. 



Proof: We follow |D2]. Let / G S and 

Ut) = e^*V-e-i|(^+-i^)/^ f^t) = e-**^*/. (3.47) 
By the fundamental Theorem of calculus, we can write 

iiUt) - fit)) = 1^ (e^^V-i^e--V-/^(,) _ K*f{s)) ds (3.48) 

Hence, 

\\fx{t)-f{t)\\ < \\K\\ f \\f,{s)-f{s)\\ ds + T{xM) (3.49) 
Jo 

12 



where 



^(x, to) = sup 
0<t<to 



(3.50) 



Now, 



SO that we can integrate ()3.5U|) by parts to obtain 



(3.51) 



(3.52) 



E 
E 
E 

3^k 



X d ^isiE,-E,)/xp,j^p^^-isK#f 



i{Ej - Ek) ds 



iiEj - Ek) 



X 



^is{E,-Ek)/xp -^p^^-isK* J 



+ 



lo {Ej - Ek) 
Hence, using ||i^*|| < ||i^||, we can bound (|3.52j) by 



E 



x\\K\\{2 + t\\K\\)e\\^\\* 
\Ej — Ek\ 



(3.53) 



Thus, 



^{x,to) < ma.x{2,to){r —r 



2 ,x(l + ||if||)||i^||ell^l!*o 



infj^fc \ Ej — Ek\ 



(3.54) 



At this point we can invoke Gronwall's Lemma, the above estimate and (|3.49|) to finish the 
proof. I 

From these two Lemmas, we immediately get the 

Proposition 3.1 Let V{x), x E [0, xq) and R be bounded operators on a Banach space B 
such that, in the operator norm, V{x) = 1^(0) + xR + 0{x'^), where V{d) is an isometry 



for any < t < to, if x ^ in such a way that t/x G 



admitting the spectral decomposition V{0) = X]j=o ^ ■'Pj and let h = 'Y^j^oEjPj . Then, 



y(0)-*/^y(x)*/^ = e*^'"^* + 0(x), in norm, 
where K* = Y.]=oPjKPj, for any K G C{B). 
Remarks: 

i) The operator in the exponent can be rewritten as 

^ihp# ^ ^ihp# ^ ^^ihp^# ^ {Re''')*. 



(3.55) 



(3.56) 
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ii) The hypotheses are made on the isometry V{0), not on the operator h. 



We can now derive our first results concerning the weak Hmit in the Schrodinger picture. 
We do so in the general setup described in HI. We further assume: 

H2: 

The restriction Hp(0) of H{0) to PB is diagonalizable and reads 

r 

Hp{0) = ^EjPj, with dim(Pj)<oo, r finite . (3.57) 

j=0 

Moreover, the operator Pe~^'^^^^^ = Pe-^'^^pi^) is an isometry on PB. 
Note that this implies Pq-'^tHp{q) -g inyertible and 

r r 

P = Y^Pj, EjEW Vj = 0, • • • , r, and Pe''^^^^^ = ^ e''^^^ Pj , (3.58) 
j=o j=o 

where the projectors Pj are eigenprojectors of Pe~''^^(°^ iff the e~''^^^ 's are distinct. In 
case i3 is a finite dimensional Hilbert space and H{0) is self adjoint, H2 is automatically 
true. 



Proposition 3.2 Let H{X) and P on B satisfy HI and H2. Further assume t > is 
such that the values are {e~*'^^j }^^Q are distinct. Then, for any <t < oo, 

lim e-*^(0)A' [Pe--^Wpl*^^'=e*r'"M on PB (3.59) 

A^O L J 

where 

r'"(r) = e'^^(°)G(T)# = - r ds r dtWe-'*^"^°^-^^^W* , (3.60) 

Jo Jo 

and K* = YJ'j PjKPj for any K G C{PB). 
Remarks: 

i) In case some values among {e^"^^Yj=o coincide, the result holds whith the P^-'s replaced 
by Ilj's, the spectral projectors of Pe^*"^^*^*^) |pg. 

ii) If, for any j = 0, . . . , r, the reduced resolvents RQ{Ej) = {H{0) — -Ej)|Qg all exist, with 
Q = I-P, then 

r"(r) = -j:P,iyPQ(£;,)(PQ(£;,)-PQ(£;,)e--(^(0)-^^)l'2«-irl)H-P,(3.61) 

i=o 

iii) If ;S is a Hilbert space, and H{X) is self-adjoint with dim Pj = 1, we can express V'" in 
yet another way. We write Pj = \^Pj){'-Pj\ and introduce d^Y {E), j = 0, • • • , r, the spectral 
measures of the vectors Wtpj = QWipj, with respect to H{0)\qb. Then, if ^ denotes the 
Fourier transform, 

r'"(r) = - V dt ^it)e''^^ Pj. (3.62) 
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Proof of Proposition 13. 2t 

As we are to work in PB, we will write Ap for PAP etc... Our assumption on r makes 
the eigenvalues of e^^'^^pi'^) distinct so that the Pjs are eigenprojectors of both Hp{0) and 
gir//p(o)_ rpj^g^^ Lemma O shows that V{x) := p(e-»^(^(o)+v^^)p, x = A^, satisfies the 
hypotheses of Proposition 13.11 with h = tHp{0), R = Gp{t) and r > fixed. Hence the 
result, making use of e~^'^^''^^^ = e~*'^^*^°^P. The last statement fohows from (|3.24j) . I 

We are now in a position to state the existence of a contraction semi-group on PB 
obtained by means of a weak limit for our specific time dependent Hamiltonian model. 
The following is a direct application of Proposition 13.21 



Corollary 3.2 Let U{n,0) be defined on Hq ^H, where Ho is separable, by \2.11[ Wl^ 

\2.(^) . let P = I and let {£'j}j=o,----r be the eigenvalues of Hq associated with 

eigenprojectors {Pj}j=o,---.r- Assume the values {6"*"^^^ }j=o,...,r cLf^ distinct. Then, for any 
fixed d <t < CO, 



lim 

A-.0 



e 



-*^W/^'PC/(t/A2,0)Pl =e*r"'M on PB (3.63) 



G{t)* = - r ds f dtY^Y. Pjy^e-''^^^+^^~^^WraPj (3.64) 

•^0 -^0 ,=Om=l 



where 

r"'(r) = e^^-^(°) 

^0 JO 

generates a contraction semi-group and # corresponds to the set of eigenprojectors {Pj }j=o,..., 
Remarks: 

0) The macroscopic time scale at which we observe the system is T = rt/}? oo. 

1) There are cases where T'^{t) generates a group of isometries. 

ii) Again, if the e~*'^^^ 's are not distinct, we have to take the spectral projectors of e~"^° 
instead of the Pj 's in the definition of the operation 

iii) Note that the effective dynamics commutes with /iq, so that no transition between the 
eigenspaces of /iq can take place. However, if the e~*'^^j's are not distinct, transitions 
between different eigenspaces of /iq corresponding to the same eigenvalue of e~"^'^ are 
possible. 

iv) In case /iq is non degenerate, r = d and we can write Pj = \xj){xj\, with Xj the 
eigenvector associated with Ej, and 







r"'W = -E EEl^^'^l^-^J-)!' / dte-^*(^*-^^+^™))|x,)(x,-|, (3.65) 

j=0 \fc=Om=l 

where the double integral equals 



Proof of Corollary ITT^ 

By Lemma l3. II above. 
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where conditions HI and H2 are met and Proposition 13.21 applies. The fact that r"'(r) 
generates a contraction semigroup in that case stems from the a priori bound, uniform in 

t,T,X, 

The expression for r'"(r) comes from the exphcit evaluation of (|3.6U|) in our model. I 



3.4 Different Time Scales 

Looking at the dependence in r of the result in Corollary 13 .21 we observe that we can ob- 
tain a different non-trivial effective evolution with our conventional weak limit approach, 
provided one further makes the time scale r ^ and, at the same time, increases the pa- 
rameter t to t/r'^. This yields a macroscopic time scale given by T = t/{T}?) oo. We'll 
come back to this point also, when we deal with the Heisenberg evolution of observables. 



Using the first expression (|3.24)) . one immediately gets 



lim lim 

t/(TA)2gN 



e*-*^W/(-")V^7(t/(TA)2,0)Pl = lime*r"'M/^' ^e*r\ 



where 



r 1 " 



(3.69) 



(3.70) 



j=0 



i=l 



Note that under the hypotheses of Corollarv 13.21 the spectral projectors of Hq and e ^"^^^ 
coincide when r ^ 0. 



This calls for a redefinition of the scaling, right from the beginning of the calculation, 
in order to arrive at the same result, without resorting to iterated limits, as above. This is 
at this point that we need to consider the dependence in r of the previous steps. 

We state below is our main theorem regarding this issue in the general Banach space 
framework under hypotheses HI and H2. Actually, the application above is a consequence 
of the theorem to come. The study at positive temperature in Heisenberg picture of the 
forthcoming Sections will rely on this result as well. 

Theorem 3.1 Suppose Hypotheses HI and H2 hold true and further assume the spec- 
tral projectors Pj, j = 0, •••,r, of e~^'^^''^^^ coincide with those of Hp{0) on PB. Set 
K* = PjKPj, for K G C{B). 

A) Then, for any < to < cxd, there exists < c < cxd such that for any < t < to, the 
following estimate holds in the limit A^r 0, A^r^ 0, and t/{\T)'^ £ N.' 



,iH(Q)t/{\^T) 



p^-iT{H{0)+\W)p 



< c(AV2 + AV). (3.71) 



B) Then, for any < to < oo, there exists < c < oo such that for any < t < to, the 



following estimate holds in the limit A r ^0, t ^ 0, and t/{\T) G N.' 



p^-ir{H{0)+\W)p 



-t{W'-)*/2 



< c{t + AV) 



(3.72) 
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Remarks: 

0) If T is small enough, the spectral projectors of e~^'^^^^^^ and Hp{0) on PB coincide. 

1) If r is fixed, part A) of the Theorem coincides with Proposition 13.21 with i := t/r"^ in 
place of t. 



Proof: We only need to consider the case r > small, where Remark 0) applies. We 
proceed in two steps, using Lemmas 13 . 41 and 13 . 51 in sequence. Let X = AV2. The expansions 
provided in Lemma 13.31 yield 

Pe-^<H{Q)+^W)p ^ ^-^rHp(0) ^ xGp{t)/t'' + 0{x^), (3.73) 
with Gp{t)/t'^ = 0(1) and reminder uniform in r ^ 0. Hence, 

||Pg-ir(^^(0)+Al4^)p|| ^ 1 + 0(3,)^ uniformly in r. (3.74) 

As e-*^-f^^(°) = e-'-'^^^Pj, with Pj independent of r, the operator J(r) defined in (l3:39|l 
reads 

•^W = E ^here a,fc(r) = ( ^"^^^^ (3.75) 



j,k=0 

Hence, 



ie J j = k. 



aikir) =i + 0(T), = -^?^ + 0(r) and J(t) = 0(1) as r ^ 0.(3.76) 

2 

Now, using ()3.21|) with coupling constant x/t (and the first remark following Lemma l3.3() . 
we can write for x/t small, uniformly in r, 

g-ir{Hp(0) + f J(r)) ^ ^3^77^ 

where the operator in the bracket above is 0(1) as r — > 0. Hence 

||e-*^(-^^(°)+r^M)|| = l + 0(a;), uniformly in r. (3.78) 
Thus, we apply Lemma |H31 to get 



p^-ir{H{0)+XW) p 



_ g-ii{rHp(0)+xJ{T)) 



0{x^), (3.79) 



as a; ^ 0, and - — > 0, with a remainder uniform in r. 



We now turn to the second step. We can write 
Therefore, by Lemma 13.51 and the last statment of 1)3. 76() . 
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uniformly in r. Hence, for any given to, we get the existence of a constant < c < oo, 
uniform in r, such that for all < t < to < c») 

<c(AV + AV2), (3.82) 

as A^T and A^r^ go to zero in such a way that t/ (Ar)^ G N, which is part A) of the Theorem. 
Part B) follows from the first statements in ()3.76|) and of the fact that the projectors P^'s 
are independent of t. I 



e 



:Hp{0) 



-itJ*{T)) 



As a direct Corollary, we get. 



Corollary 3.3 Let U{n, 0) be defined on Tio (8) fi, Tio i separable Hilbert space, by \2.1% 
\2.1(A rOI) . let P = 1® and let {£'j}j=o, the eigenvalues of ho associated with 

eigenprojectors {Pj}j=o,---,r- Then, for any <t <to, 



lim 

t-.0,a2t-.0 
t/(TA)26N 



^*^(0)/(^^)Vc/(V(rA)2,0)p' 



(3.83) 



where Tq* = T!j=oPj^oPj, and Tq = -i ELi ^^^Z- 



4 Heisenberg representation for non-zero temper- 
ature 

From now on, we stick to our model Hamiltonian characterized by hypothesis HO. We first 
express the evolution at positive temperature of observables B of the small system ()4.4() 
after k repeated interactions as the action of the /c-th power of an operator Up{X, r) on TLq. 
This reflects the Markovian nature of our model. 

This is done in Proposition 14. II This allows us to apply Theorem 13. II again to compute 
the weak limit in Theorem 14.11 Let us mention here already that we perform a complete 
analysis of the special case where both the small system and the individual spins of the 
chain live in in the last Section of the paper. 

Let us deflne the equilibrium state lo{I3)n of a chain of N spins at inverse temperature 
/? by a tensor product of individual diagonal density matrices of the form 



1 + E ■=! e 



-135, 



(I 




Vo 



e 



\ 



,—f3Sa*a 



(4.1) 



in the basis {lo, xi, • • • , x„} of C" , i.e. 

'^(/3)iv = r{(3) r{/3) O • • • r(/3). 



(4.2) 



The individual density matrices r(/3) are deflned by Gibbs prescription for the Hamil- 
tonians at each site 



(4.3) 



i=l 
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corresponding to our model H2.7|) 

Our spin chain is of finite length N, but, as we will see below, only the first k spins matter 
to study the time evolution up to time k. This will allow us to take the thermodynamical 
limit by hand. If p is any state on C^"*"^, the initial state of the small system plus spin 
chain is p uj{[3)n- We shall study the Heisenberg evolution of observables of the form 
B0Ih, where B e Md+i(C), defined by 

Bfs{k,X,T) = TTn{{I(^u;Nm U{k,0)-\B ^In)U{k,0)), (4.4) 
where, for any A £ OJ-Lq Ti), 

Tr-^(yl) = ( ^^(xi ® xs\Axj ® xs) I with xq = uj, (4.5) 

denotes the partial trace taken on the spin variables only. Hence, the expectation in the 
state p of the observable B after k interactions over a time interval of length kr with the 
chain at inverse temperature (5 is given by 

{B{k, f3))p = TVc.+i {pBp{k, A, r)). (4.6) 

Remark: 

In case TLq is infinite dimensional, the definitions (|4.4|) and (|4.5jl hold, mutatis mutandis. 
For instance, consider B G C{Ti.o) in H4.4() . where (|4.5|) should be read as 

Tvn{A) = Y,{^no(^xs\AIno(S)\xs), (4.7) 
s 

with a slight abuse of notations. 

4.1 Markov Properties 

Recall that 

U{k, 0)-\B ® In)U{k, 0) = C/i*C/| ■■■U;{B0 InMUk-i ---Ui, (4.8) 

where Uj is non-trivial on C"'^"'^ C""*"^ only. 

Let us specify a bit more the partial trace operator Tv-)^{{I u)n{i3)) A), where A is an 
operator on C^+i ® Uf^^C]+^. 

Lemma 4.1 Let us denote the matrix elements of A as follows 

^S,S' = {xi(E>Xs\Axj(g>Xs'), (4.9) 

where i,j belong to {0, • • • , d}, and S, S' run over subsets of {{1, • • • , N} x {1, • • • , n}}^ 
as in E^) Then 

^ e-I^T:?=iii\s\i 

Trniil^ivM) ^kj = T.j^^^^n—z^A's's (4-10) 

s 

where, for 

S = {(A:i, ii), {k2,i2), • • • , {km,im)} C (N X {1, 2, • • • , n}r (4.11) 
with all 1 < kj < N distinct and m = 0, ■ ■ ■ , N , 

\S\l = #{kr S.t. ir = l}. (4.12) 
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Proof: Follows directly from 

(1 + E,- e-^^^)^^^"" " (1 + E,- e-^^^) 



We now further compute the action of U{k,0) given by the product of UjS. Let us 
denote the vectors u (8) and xj Xg by no ® n2, • • • , tin) = no |n), where no € 
{0, 1, • • • d}, and n^ G {0, 1, • • • n}, for any j = 1, - ■ ■ N, with a; c± and Xk — k and 

^{(l,ni),-,(Af,n^)} =^ 1^)- 

Recall that 



■'J 



(4.14) 



where e ^"^^^ is diagonal. More precisely, with the convention So = 0, 

e ■'no <8) |ni,n2, • • • jJ^iv) = e no <8) |ni, n2, • ' ' > J^w)- (4-15) 



Lemma 4.2 Denoting the k -independent matrix elements ofe ^'^^''{i^d+i^^ri+i = Uk\i^d+i^cn+i 
by 

d' = {nC^mlUkn' m'), (4. 16) 

we have for any N > k 

UkUk-i ■ ■ ■ U2U1 no |ni, • • • , tin) = (4-17) 

w;/tere 

j=l \j<l<N Kj J 

Proof: Consequence of the iteration of formulae of the type 

Ui no ® |ni, • • • , njv) = ^ ^"^ f^TO ® |mi, n2, ng, • • • , niv).(4.19) 



A consequence of these formulae is that we can consider spin chains consisting in k 
spins only: 

Lemma 4.3 For any N > k, 

Trn{^ ^ ^NiP) UlU^ ■■■Ul{B® ln)UkUk-i • • • C/i) = 

Trnil ® uJkiP) U^UI ■ ■ ■ U^B In)UkUk-i ■ ■ ■ C/i)(4.20) 
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Proof: Obvious from the tensor product structure of uJ]\f{P). 
To proceed, let us adopt the fohowing block notation 



U 



-iT{H{0)+\W) 



where 



U„ 



' m,m' 



0,1 



m,m' 
1,1 



0,d 



m,m' ^ 
l,d 



m,m' m,m' 



u 



In terms of the notations of the previous Section, 

(PUP PUQ\ 
\QUP QUQj' 

we have the identifications 

PUP^Uo,o, QUQ^ : •.. : 

\f^i,l ••• Un,n J 

PC/Q~(C/o,i ••• C/o,n), QUPc^iUifl 



Let us finally denote the inverse of f7 = (U^^,) by 

y = iKi') = = {U-'Z^,) G M(i+,)(i+„)(C), 
so that we have for any m and n 

Km = ym,n G Mi+d(C). 

With these notations, we have 



Un,o , 



(4.21) 



(4.22) 



(4.23) 



(4.24) 

(4.25) 
(4.26) 



Lemma 4.4 The matrix elements ofU{k,0) ^ {B^I-)^) U{k,0) in the orthonormal basis 
{no |ni, • • • , Uk)} = {riQ (gi |n)} read 



(no O n| (Uk ■■■Ui)*B(^ In (Uk ■ ■ ■ Ui) hq (g) n) 



(4.27) 



-iT(if>{0,n)-ip{0,n)) 



■ Uqjii ,ni ) 



no, no 



fne{0,---,n}* 



Proof: Expand the products and make use of Lemma 14.21 and (|4.18|) . 

The above Lemmas and (|4.4() lead us to the study of the matrix in M^^i (C) 



E 



m — (m-j^ , ■ ■ - m^, ) 



(4.28) 
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in various limiting cases as A and/or r go to zero, with the notation 

\n\i = f,{nr s.t. nr = l} = \S\i. (4.29) 

We introduce operators on the Hilbert space M^+i (C) equipped with the scalar product 
{A\B) = Tr{A*B), for any A,B e Md+i{C) by 

Um,m'{A) := Vm'^mAUm,m', (m, m') G {0, 1, • • • , n}^. (4.30) 

These operators are linear and one has with respect to the above scalar product, 

The composition of such operators will be denoted as follows 

^^m' ,n' ^m,n{A) = Vn' ,m'^n,m AUm,nUm' ,n' ■ (4.32) 

We are now in a position to express the Markovian nature of the evolution of our 
observables: 

Proposition 4.1 In terms of the operators defined above, we can write 

+Uifl + e~^^^Ui,i + ■■■ + e-^^-Ui^n 

\ k 

-(^^17 J , _L . . . _1_ c-/3'5n7 



^Up{\Tf{B). (4.33) 

Proof: By definition of hlm,n we have 



Furthermore introducing 3^m,n = e~^"f^Um,n, we get 

Bp{k,X,T)= Yl y^,,n,---y„,{B). (4.35) 

m = {mi,.--,mj.) 

There are (n + 1)^ distinct operators ym,m' in that expression, and the set of vectors n, m 
in the sum yields all different ways of composing k of them. Therefore 

Bp{k, A, r) = (Jo,o + • • • + Jo,n + • • • + 3^n,0 + • • • + yn,nf {B).{A.36) 



Remark: 

The formula of Proposition 14.11 holds if TCq is a separable Hilbert space, provided the 
decomposition of operators A in (|4.21j) is interpreted as Apg £ CiTio), q,p £ {l,---,n}, 
with 

Apg = Iho ® \p){p\AIn, \q){q\, (4.37) 
and the identification Tio <Si C\q) ~ Tio, for all q. 
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4.2 Weak Limit in the Heisenberg Picture 

The A-dependence in Bf^{k, X,t) comes from the definition 

[/ = C/,(A) = e-^-W0)W)^ (438) 

which imphes that the lAn^mS depend on A as well, in an analytic fashion, and will be de- 
noted Un,m{X)- Expliciting the A dependence in Bp{k^ A, r), the weak limit corresponds to 
taking k = t/X"^ and computing the behavior of Bp{t/X^ , A, r), as A — > (keeping r fixed). 
We shall use the same strategy as in the previous Section and Lemma to identify the 
weak limit by means of perturbation theory. We shall also eventually consider the possi- 
bility of letting r — > 0, therefore we explicit the behavior in r of the expansions below. 

Consequently, with (|4.24|) and Corollarv 13.11 we get 



Lemma 4.5 Let U be given by \4.^^ , with H{0), W self adjoint and satisfying HI, and 
further assume H{0) is diagonal with respect to the basis it^.5|) . If tlm,m'{^) is defined by 
i4-3^ As X ^ 0, we get the expansions 



UoflW = UoM + xHfi + 0(AV4) (4.39) 

llm,m'{X) = ll^,m'{0) + X^uj^^m' + OiX^T^), m,m' > 1 (4.40) 

Uo,m{X)=X^U^']l + 0{X\''), m>l (4.41) 

Z^„,o(A) = A2z^i'i + 0(AV4), m>l (4.42) 
where, for all < m,m' < n 

Um,m'iO)iB) = <5™,™.e^"^--(°) B e— (4.43) 

= W'(G„^,™(-T)Se-^-^™-W + e^^^--(°)5G™,^(T)), (4.44) 

and, for all 1 < m, 

^oil(^) = Frn,o{-T)BFo,m{r), (4.45) 

uj^liB) = Fo,^(-r)i?F„,o(r). (4.46) 
This Lemma allows us to perform the analysis of the operator defined in Proposition 

Up{X,T) = Z{py^ Ui^mWe-^"-" , as A ^ 0, (4.47) 

0<m<n 
0<;<n 

with the convention = and Z{f3) = Sj=o ■ Recall that 

Bp{k,X,T)=Up{X,Tf{B). (4.48) 
Moreover, using the fact, see (|3.4j) . 

F^,™(0) =Fo,o(0)+5„^~/^o + <5m, (4.49) 
we get for all < m < n 

Um,n.mB) =Uo,omB) ^ e'-f^^Be-'^^'^ = e'^^^^^'\B). (4.50) 
We have thus shown the 
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Lemma 4.6 Assume the hypotheses of Lemma \4-5\ Then 

A2 



+ 0(AV 



_m=l 

= Z^o,o(0) + A2^(/5)-^r^ + 0(AV), (4.51) 
with Tf3 = T(s{T) = 0(t2). 

The above operator enjoys the following symmetry property 
Lemma 4.7 For any B £ Md+i(C), 



Tr{BTp{B*)) = Tr{B*Tp{B)). (4.52) 
Proof: Due to 

Fn,m{-T) = Fm,n(T)*, m^U (4.53) 

G„,„(-r) = G„,„(r)* (4.54) 
and to the structure of T/j, the result will be proven once we show that for all A,B,C £ 



Tt{B*ABC + B*C*BA*) = Tt{BAB*C + BC*B*A*). (4.55) 

But this follows from Tri? = Tri?-'", where denotes the transpose, and from the cyclicity 
of the trace again. I 

Recall also the property 

Up{X,T){I)=I r^(I)=0, (4.56) 

and the fact that in case the spectrum {£'j}j=o,---,d is non-degenerate and {|xj)}j=o,---,rf 
denotes the corresponding eigenvectors, the unitary h(ofi{0) has degenerate spectrum: 

Uo,om\xj){xk\) = e^^(^^-^'=) \xj){xk\, yO<j,k<d. (4.57) 

That is, (t(Z//o,o(0)) = {e^'^^^^~~^''^}o<j,k<d, so that 1 is (i+ 1 times degenerate at least. 

We are in the same position as in the proof of Proposition H3.2() . Therefore, we can 
compute the weak limit from Proposition 13.11 immediately to get the following 

Theorem 4.1 LetUp{\,T) he given by UT^ , andUofi{0), by |737p . Lei {e*^^'}i=i,...,,. 
he the set of distinct eigenvalues ofUofi{0) and denote by Pi the corresponding orthogonal 
projectors. Then 

Ikn Uo,o{Or'/^'Bf3{t/X^,X,T)= (4.58) 
Ihn Uo,o{Or'^^'m^,rf/^\B) = e'^'^{B), 

t/A2eN 



were 



r;^(i?) = ^ {1^0,0(0)-' Tp)* (5), (4.59) 
with # corresponding to the set of projectors 
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Remarks: 

0) In order to make the generator completely explicit, one needs to analyse the prop- 
erties T^, i.e. of the operators Vj defining the coupling, within the eigenspaces of Z^o,o(0)- 
A non trivial example is worked out in Section 6, see Proposition 16.11 

1) The degeneracy of the eigenvalue 1 of Wo,o(0) is responsible for the existence of a non- 
trivial invariant sub-algebra of observables which is the commutant of /iq. 

As in Section 2, we generalize our result to the regime A^r — > 0, r ^ 0, by switching 
to the macroscopic time scale T = t/{\^T) — > oo. We first compute 

^^^^^ = }^ ^''zZ)r^^^ ^^^ ^ "2^^^'°'°^ + + ^^-^^^ 

ra=l ^ ^ ' ^ 

which, using the following formulas for m > 1 

n 

to express the operators Wmm' in terms of Vm, eventually becomes 
1 " / 1 

TM^) = ^^T.e-^''-{ymBv^--iv^v:,B+Bv^v;:,: 



Zip) ^ > 'm 2' 

^ ' m=l ^ 

^V*^BYm - \{V^V^B + BV;,Vm). (4.62) 

We note here that this operator has the form of the dissipative part of a Lindblad generator. 
We'll come back to this operator Tp in connection to the modelization in terms of Quantum 
Noises proposed in |APj and |LMj . in the next Section. 

Corollary 4.1 Assume the hypotheses of Theorem \4.1\ Then with t/{TX)^ = A; G N, 

lim ^/o,o(0)-*/(^^)'i?;3(t/(TA)2,A,r)) = (4.63) 

t/(TA)2gN 

lim Z^o,o(0)-*/(^^)'z^^(A,T)*/^'(i3) = e'^^*{B), 

t-»0,a2t^0 
t/(TA)2gN 

were Tp{B) is defined in \4.6l3{j . 

Proof: We can simply repeat the arguments of the proof Theorem 13.11 once we note the 
following facts: i) The operator Z^o,o(0) = e*'^'''"''] is unitary on Mrf+i(C), with spectral 
projectors that are independent of r as r — > and eigenvalues of the form e*"^^^ . ii) 
Introducing x = (Ar)^, (|4.51|) states that uniformly in r, 

Up{\T) =U^fl{Q) + xTp{T)/{T^Z{())) + 0{x^), (4.64) 

where Tp{t)/t'^ — > as t — > 0. I 
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4.3 Evolution of states 

Let us close this Section by briefly recalling some consequences of these results about the 
evolution of states, i.e. trace one positive matrices. This is conveniently done in our setup 
by using duality with respect to the scalar product {A\B) = Tr{A*B). 

If r is the generator of the dynamics of observables, B is an observable and p is a state, 
then for any t G M, 

Tr(pe*^(S)) = Tr(e*^* {p)B) (4.65) 

where the generator of the dynamics of the states is T^, such that for all states p and 
observables B, 

TT{{r,{p)rB) = {p\r{B)) = {r*p\B). (4.66) 



In the particular case where the observables Pjk = \xj){xk\, with the notations of 
form an orthonormal basis of eigenvectors of the restricted uncoupled evolution Uq^q, the 
corresponding eigenprojectors are denoted by Hjk and act as 

Ujk{B) = PjkTr{\xk){xj\B) = P,k{xj\Bxk)no, (4-67) 

where the subscript TCq denotes the scalar product within TCq. Hence, the # operation on 
the operator T with respect to the projectors Jljk is given by 

T*{B) = Y,^jkrUjk{B) = \xj){xk\ {xj\r{\xj){xk\)xk)no{xj\Bxk)no- (4-68) 

j,k j,k 

Therefore, one computes that the corresponding generator of states, {T"^)* is given by 

rf (p) = ^ \xj){xk\ {xk\r{\xj){xk\)*Xj)no{xj\pXk)no- (4-69) 

j,k 

Consequently, 

(r#), = Y,^jkr,Ujk = (r.)#. (4.70) 

j,k 

We note that states defined as functions of the Hamiltonian ho of the small system 
form an invariant subspace of sets whose Markovian dynamics is characterized by the 
scalars {{xj\r{\xj){xj\)xj)Ho}j=o,-,d- 



5 Beyond the perturbative regime: A^r = 1 

We consider here the regime A^r = 1, and r — > used in |APj in their construction of the 
field of quantum noises. It can be viewed as a regime where the weak limit scaling holds 
at the microscopic level, while, at the macroscopic level, T = ^/(tA^) is kept finite. 

As we saw in Corollaries 13.31 and 14.11 in the Schrodinger and Heisenberg pictures re- 
spectively, the small parameter that allows to make use of perturbation theory to compute 
the effective evolution is the combination A^r. Therefore, we have to resort to a differ- 
ent technique since our scaling imposes a non-perturbative regime. Our main tool will be 
Chernoff's Theorem as we now explain. 
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5.1 Schrodinger Evolution 

Let us start with the Schrodinger effective evolution under the following assumptions: 



H3: Hypothesis HI holds with B a Hilbert space and P, H{X) = H{0)+XW self-adjoint. 

In the scaling adopted here, the number of interactions n has to grow like n = t/r. This 
is in keeping with by the fact that in all cases considered so far, n = t/{\T)'^ = t/r. Note 
that the macroscopic time T = rn = t is finite here. Therefore, according to the analysis 
of Section 3, we are led to study 

PU{t/T,0)P = rpe-*(^^(°)+v^^)p]*^^ , as T ^ 0, t/r G N*. (5.1) 

This limit is easily computed by applying the following version of Chernoff 's Theorem, see 
e.g. jBRj . jD2] or [Paz , which suffices for our purpose: 

Theorem 5.1 Let S{t) defined on a Banach space B be such that S{0) = I, and \\S{t)\\ < 
1, for all r > 0. //, limT-_^o ■'"^^(5'(r) — I) = F in the strong sense exists in C{B) and 
generates a contraction semi-group, then 

s-lim5(t/n)" = e*^. (5.2) 

Now, it is easily checked that 

5(r) := Pe-^(^^(°)+^^)p on the subspace PB (5.3) 

satisfies the first requirements. Then, by expanding the exponent and making use of the 
properties of -ff (0) and W , we can write 

S{t) = {l- itH{Q)p - ^-{W^)p + 0(r2)) . (5.4) 



It thus implies 



S'{T)\r=o = -iH{<d)p - = r G C{PB). (5.5) 

Now r is dissipative, since Mip G PB 

^{<f\T^) = -^{ip\PWQWP^)/2 = -\\QWP^\\b/2 < 0. (5.6) 

Hence, by Lumer-Phillips, see |Pazj . F generates a contraction semigroup. Therefore 
Theorem 5.2 Under the hypothesis H3, for any t > fixed, 

s- hm PU{t/T,0)P = s- lim [pe-(-^(o)+v^^)pl = e-*(*^W^+^). (5.7) 

T—>0 T—>0 |_ J 

t/rGM t/rGM 

Remark: Specializing to our model Hamiltonian, we get that the effective dynamics on 
PB is 



(5.8) 



Apart from the self-adjoint part Hq stemming from the uncoupled evolution, the main 
difference with respect to the corresponding weak coupling result in Corollary 13.31 lies in 
the absence of the # operation on the dissipative part ^ J2j ^j'^j of generator. This 
prevents the spectral subspaces of Kq from being invariant under the effective dynamics. 
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5.2 Heisenberg Evolution 

Let us now turn to the more interesting case of the Heisenberg dynamics of observables when 
the spins are at equihbrium at inverse temperature (3. We assume the general hypothesis 
HO, i.e. we stick to our matrix model, even though certain results below hold for more 
general situations. 

The analysis of Section 4 shows that the evolution of an observable B G M^_(_i(C) after 
k repeated interactions reads 

B^Bp{k,\,T)=Up{\Tf{B) (5.9) 

with Upi^X, t) defined by (|4.47|) . where we explicited the dependence in r in the notation. We 
want to apply Chernoff's Theorem again to the operator valued function r i— > Up{l/ y/T,T) 
on £(Afrf_|_i(C)). In order to check the first hypotheses we recall the formula (see ()4.1U() 1 

Up{X,T){B) = Trn{{l(^iom)U-Hl,O){B0l)U{l,O)) (5.10) 

where B(t),, = ([/-^(l, 0)(S I)C/(1, 0))q5 = ^^[/-^(l, 0)(B (g) I)C/(1, 0)Pg according to 
the block notation 1)4. 21() . with the corresponding orthogonal projectors Pq. Identifying 
p^(^{n+i)(d+i) ^-^j^ ^ (j-d+i^ deduce from the above formula that Z^/3(A,t) is a 
contraction for any value of the parameters: 

n 

\\Ufs{X,T){B)\\no < E Wll^(^)'?'?ll^o (5-11) 



q=0 



q=0 

n 



q=0 



Moreover, Uf3{l / y/r , T)\r=o = I, so we are left with the computation of the derivative w.r.t. 
T at the origin. This involves the control of the operator Ur{X) 1)3.20(1 as r — > and 
A = 1/\/t — > oo, as in the previous paragraph. Let us get estimates in a more systematic 
way than above. So far, all our estimates are derived for both A and r going to zero or at 
most finite. However, the expansion of Ur{\) in powers of A is convergent, with r dependent 
coefficients we control sufficiently well. Indeed, ()H. ;-};-{ |) yields 

Ur{X) = e-*^^We(A,T) = Ee-^"''^°^e„(A,r), (5.12) 

n>0 

where 0„ contains n operators W and satisfies 

||G„(A,r)|| =0((rA)7n!). (5.13) 
Using the fact that (Ar)" = r*^/^ and that W is off-diagonal with respect to P and Q, 



we get that the replacement of A by l/yr doesn't spoil the estimates as r ^ given in 
Proposition 14. II and Lemma 14.51 Those together with the computation (|4.62() yield 

Up{l/V^,T){B) = e'-^''Be-'-''^' + {Z{P)T)~%{T){B) + 0{T^) 



e 



irhoQ^-irho ^ rVpiB) + 0(t2), (5.14) 
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where, see 1)4. 62() . 

m=l ^ ^ 

- l{V:,VmB + SK;^.). (5.15) 

Hence, the derivative at the origin exists and is given by 

U(,{l/V^, r)'(i?)|,=o = i[ho, B] + Tf,{B). (5.16) 

We recognize at once that Tj3{B) is the dissipative part of a Lindblad operator of the 
form 

2m 

E ^^^^J - 2 ^^^^^^ + ^^^-^i) (^-^^^ 

i=i 

with 

L,- = Vj, l<j<m and L,- = V* , m + l< j < 2m. (5.18) 

By the Theorem of Lindblad, see e.g. |AFj . we know that 

i[ho,B]+Tp{B) (5.19) 

generates a completely positive semigroup of contractions. Therefore, we are in a position 
to apply Chernoff 's theorem to eventually get 

Theorem 5.3 Assume hypothesis HO where TCq is a separable Hilhert space and ho, the 
Vj 's and B are bounded on TLq. Let Bpit/r, 1/\/t, t) be defined by i4.4\ ), ^/^(A, r) is defined 
by proposition \4-.l\ and the Remark following it. Then 

s - lim Bs{t/T,l/^/T,T) =s - lim ^/«(1/V^, r)*/"(5) = e*(^['*«'-l+r/3(-))(5) (5.20) 

T — tO T — 

t/reN t/T6N 

with a Lindblad generator i[ho, ■] + r^(-) explicited in i5.17\ ) 
Remarks: 

i) Let us make a comparison of the above with the results of |APj . Section IV. 2, which 
concern similar generators as ours. More precisely, (|2.6|1 corresponds to a particular case of 
the Hamiltonian of eq. (15) in |APj . with Dij = 0, Vi,j. In (AP,, the choice of time scale r 
and coupling A is such that A^r = 1, r ^ 0. A supplementary structure is present in that 
work which consists in making the suitably renormalized spins forming the chain merge 
in the limit r ^ to yield a heat bath represented by a Fock space of quantum noises. 
The limit r — > performed in the language adopted in |APj exists and yields a quantum 
Langevin equation for the whole limiting system consisting in the original small system in 
interaction with a field of quantum noises. When restricted to TCq, the effective dynamics 
of observables at zero temperature corresponds to a contraction semigroup generated by 

roo(-) = i[ho, ■] + E • - \{V:,Vm ■ + ■ V:,V^)] , (5.21) 

which coincides with Theorem 15. 31 at (j = 00. 

ii) The generator T^^ coincides with the generator 1)4. 62() obtained in Corollarv 14.11 in the 
scaling A^r — > 0, r — > 0, modulo the # operation, which appears as a trade mark of the 
perturbative regime. 
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5.3 The Continuous Limit 



For completeness, we mention here the easier cases of continuous Umit characterized by 
r ^ and A constant. The omitted proof are quite analogous to those of the previous 
Section. 

First considering the Schrodinger picture, we get 
Proposition 5.1 Assume the hypothesis H3 holds and fix A = 1. Then, 

s- lim PU{t/T,0)P = s- lim \pe-'<^^^^+^^ pY^^ = e''^^'^^^'' . (5.22) 

T— *0 T — '0 |_ J 

t/T6N t/reN 

The Heisenberg evolution also yields a unitary effective evolution in the continuous limit: 
Proposition 5.2 Consider the matrix model of Section 2 and fix A = 1. Then, 

lim BM/t,t, 1) = e^[^0'-l(B). (5.23) 

T — >0 

t/rGN 

In order to make explicit the results of Sectional we provide below a detailed analysis 
of the case d = n = 1. 

6 The case d = n = 1 

In that Section, we focus on the first non-trivial case where the small system lives on 
and the heat bath is formed by a chain of spins 1/2. We provide explicit formulas for 
Tfj and which are valid for any coupling operator V appearing in ()2.8|) . We further 
diagonalize the restriction of T^ to the degenerate subspace Ker(Z^o,o(0) — 1) in order to 
determine the subalgebra of observables invariant under the effective dynamics in the weak 
coupling limit (keeping r fixed). 

For Tio = 'C?,n = ®j>i<C?, we write for t G [r/c - 1, Tk[ in Hq ® C^, 

H{t,X) = H{\) = H{0) + \W, (6.1) 

where 

H{0) = ho 01 + 10 6a*a, W = V*0a + V0a*. (6.2) 

We choose, without loss of generality, /iq = ecr^, e 7^ 0, so that we have in the ordered basis 

{w (X" tj, X w, (8) X, X (8) x} 

H{X) = (^^^Y gj^^^^ ^ ' with the convention ~ ^ ■ (6-3) 

Specifying the results of the previous sections to the case under study, we can write, 
uniformly in (3, as A ^ 0, 
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with 

Tp{B) = Fo,i{-T)BF,,o{r) + Go,o(.-T)Be~'^''oMo) + ^^rHoMo)BGo,o{T) (6.5) 

We use the norm induced by the scalar product {A,B) = Tt{A*B), i.e. the Hilbert- 
Schmidt norm. As easily verified, an orthonormal basis of eigenvectors for the unitary 
operator Z^o,o(0)(-) = e*'^'^'^^ • e"*'^'^'^^, with associated eigenvalues, is provided by 

{i,a„cT_,a+} ^ {l,l,e-2-^e2-^}, (6.6) 

where 



(6.7) 



^+={1 o)> ^- = (o J)' S = VV2 and a, = ajV2. 
Let us compute restricted to the subspace Ker (Z//o,o(0) — 1) appearing in T^. 
Lemma 6.1 With respect to the orthonormal basis {I, a^}, and with the notation 

we have 
w/iere 

Tp,a = -(l|i^iTll' + ^-'"ll</'ll') = -|l^iTf(l + <0- (6-10) 

a b 



Furthermore, if in i2.8\) V ^ ^ 

/ e'^^'-^^ e'^^dsa e'^^'-^^ e^'^^-'^'Ush\ 

Proof: The first column is proportional to Tp(I) = 0. The second column of the matrix is 
given by 

1 / Tr(r,(a.)) \ 

where, dropping the positive argument r in F and further making use of (|3.,36|1 and (|3.37|) . 

+ e"'" (p^.a^Fo^i + G^a,e-*^^i.^(0) + e*^^M(o)^^Gi,i) . (6.13) 



Further making use of the cyclicity of the trace, [fx^, ff„^„(0)] = 0, o"^ = I and of (|3.36|) and 
()3.37() again, we can write 

TT{a,Tp{a,))) = Tr((7,Fi*o^.^i,o) - TV(Fi*oFi,o) 

+ e-'^^(Tr(a,Fo*ia,Fo,i)-Tr(Fo*iFo,i)). (6.14) 



31 



Explicit computations on 2 x 2 matrices yields the first equality in H6.1U() . Let us turn to 
(|6rTT|) . From the definitions (|3?23|) and (jT8|) . we have 

^^"^=Ui,oW o J' (^-^^^ 

where 

Fo,i(r) = -i r e-ii^-^)Ho,o{o)v*e-''^^-^''^Us, (6.16) 

JO 

Fi,o(r) = -i r e-i(^-')Hi,iio)ve-''^o,omds. (6.17) 



By explicit computations with V as in the statement, we obtain 

which yields the expression for Tj^. . . 
By similar manipulations we get 

TT{ITp{a,)) = TY(Fi>,Fi,o) - TV(a,Fi*o^i,o) 

+ e-^^(Tr(Fo*ia,Fo,i)-TV(a,Fo*iFo,i)). (6.20) 
Now, for any F E M2(C), 

Tr(a,(FF* - F*F)) = 2{\F2i\^ - \Fu\^), (6.21) 
so that we get the first hne of 1)6. 10|) . I 

We also need to compute Tr(cj_T^((T+)) and Tr(cr_|_T^((T_)) to get T^. 
Lemma 6.2 By explicit computation and Lemma \4- T[ we have 



Tr{a.Tp{a+)) = Tr{a+Tp{a.)) 



= (i^i,o)i,i(i^i,o)2,2 + e*"^((Go,o)i,i + (^0,0)2,2) (6.22) 
+ e-'^ ((Fo,i)i,i(Fo,i)2,2 + e'^'{e'^'{Gi,i)i,i + e-^^'^(Gi,i)2,2)) • 



It remains to diagonalize the restriction of Tp to span(I, a^) to have a complete descrip- 
tion of the generator of the effective evolution. Introducing 

= ^/3l,2' ^ = ^/?2,2' {6.23) 
we actually get by perturbation theory. 

Lemma 6.3 Assume er ^ Zvr. Then, for A > 0, there exists a continuous set of eigenpro- 
jectors and eigenvalues ofhlp{X,T) denoted respectively by {Iij{X)} j=i^...^i and {uj(A)}j=i^...^4 
such that 

ni(A) = l + 0(A^), 

U2(A) = l-A2||FO^f + 0(A^), 

u^{\) = e^-^ + \^Tr{a^Tp{a+)) + 0{X^), 

Ui{X) = e-^'^' + \^Tr{a+T^{(j^)) + O(A^), (6.24) 



32 



and 

U3{X){B) = Tr{a^B)a+ + 0{\^), Tli{\){B) = Tr{a+B)a^ + 0{X^). (6.25) 

Moreover, Hq := ni(0) +112(0), 113(0) and 114(0) are the spectral projectors o/Wo,o(0) and 
{nj(0)}j=i^...^4 are those ofTp. 

Hence, we obtain the 

Proposition 6.1 Let t/}? = A; G N, and consider the Hamiltonian Then 

limZ^o,o(0)-*/^'i?(t/A2,/3, A) = lim Z^o,o(0)-*/^'^/MA, r)*/^' (i?) = e'^^B), (6.26) 

A — *0 A — ^0 

were 

+ e-2-^Tr(a_r^(a+))n3(0) + e^'^'Tria+Tp{a^))U,m. 
The dynamics of any observable is thus fully determined from these formulas. 
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